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Perturbation theory for dark matter clustering has received a lot of attention in recent years, but 
its convergence properties remain poorly justified and there is no successful model that works both 
for correlation function and for power spectrum. Here we present Halo Zeldovich approach combined 
with perturbation theory (HZPT), in which we use standard perturbation theory at one loop order 
(SPT) at very low k, and connect it to a version of the halo model, for which we adopt the Zeldovich 
approximation plus a Pade expansion of a compensated one halo term. This low k matching allows 
us to determine the one halo term amplitude and redshift evolution, both of which are in an excellent 
agreement with simulations, and approximately agree with the expected value from the halo model. 

Our Pade expansion approach of the one halo term added to Zeldovich approximation identifies a 
typical halo scale averaged over the halo mass function, the halo radius scale of order of IMpc/h, 
and a much larger halo mass compensation scale, which can be determined from SPT. The model 
gives better than one percent accurate predictions for the correlation function above 5Mpc/h at all 
redshifts, without any free parameters. With three fitted Pade expansion coefficients the agreement 


in power spectrum is good to a percent up to fc ~ 
by adding higher order terms in Pade expansion. 

PACS numbers: 98.80 

I. INTRODUCTION 

One of the major puzzles in cosmology is understand¬ 
ing the nonlinear formation of structure in the universe. 
The current state of the art are the N-body simulations, 
which have been verified to give reliable answers at 1% 
level in the power spectrum up to fc ^ Ih/Mpc [1], but are 
expensive to run, require large allocations, and are often 
not fully convergent using typical current generation box 
size and resolution. The convergence properties become 
a lot more difficult to achieve for higher order correla¬ 
tions. For example, recent studies have shown that to 
reach one percent convergence on the covariance matrix 
one needs to simulate a volume in excess of lOOOGpc^ [2], 
and the covariance matrix depends both on the cosmo¬ 
logical model and on the size of the survey one is simulat¬ 
ing, making the numerical solution to the full problem an 
excessively demanding task given the current typical re¬ 
sources. Even more importantly, we want to use cluster¬ 
ing statistics to extract information about our universe, 
and simulations do not provide much insight to the ques¬ 
tions such as where is the information content and how 
to optimally extract it from the data. 

Alternative approach is to use perturbation theory 
(PT), of which the two most prominent examples are 
standard PT (SPT) and Lagrangian PT (LPT) (see [3] 
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Ih/Mpc, which can be improved to arbitrary k 


for a review). In this approach one assumes density per¬ 
turbation is less than unity, and one expands the nonlin¬ 
ear equations perturbatively. Since density perturbations 
grow with wavevector k, this approach breaks down for 
some k > fc„i. The current approaches typically under¬ 
predict power in LPT and overpredict in one loop SPT. 
For example, we know that Zeldovich approximation is 
quite successful in modeling the correlation function on 
large scales, including the baryonic acoustic oscillations 
(BAO) wiggles [1]. However, it fails miserably in the 
power spectrum, underestimating the power at all but the 
lowest values of k, and giving predictions that are even 
below the linear theory at 0 = 0 [S]. Various one loop 
LPT extensions (one loop LPT [5], CLPT [7], CLPTs 
m somewhat improve this for the power spectrum, but 
make things worse for the correlation function [5]. Phys¬ 
ically, the problem of Zeldovich approximation and its 
extensions is that while it can describe properly the ini¬ 
tial streaming of dark matter particles, it fails to account 
for their capture into the dark matter halos: instead, the 
particles continue to stream along their trajectories, set 
by initial velocity in Zeldovich approximation, leading to 
an excessive smoothing of the power. 

Standard Eulerian PT (SPT) has the opposite prob¬ 
lem. In power spectrum it quickly overpredicts the 
amount of power at one loop level, specially at low red- 
shifts. The reason for this is that the loop integrals ex¬ 
tend over all modes, including those that are in the non¬ 
linear regime k > fcni- These modes are not in PT regime 
and typically these contrubutions are strongly suppressed 
in the simulations relative to PT, leading to too much 
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power in SPT relative to simulations. Its Fourier trans¬ 
form, using a gaussian smoothing to obtain a convergent 
correlation function, results in a worse model than the 
Zeldovich approximation around BAO, but is otherwise 
comparable. Recent work emphasized this point in the 
context of effective field theory (EFT) [8]. For example, 
for power law power spectra these integrals can be di¬ 
vergent, so PT is wrong in such situations [3]. Instead, 
it is argued that the best that one can do is to intro¬ 
duce effective field theory parameters that describe the 
correction from the small scale physics. For k < 1/Rh 
these terms can be expressed as a low k limit of PT, but 
with free parameters. A priori it is unclear what is Rh, 
and how large these corrections are for our universe. In 
particular, realistic CDM power spectra have the shape 
where most of the low k limit loop integral comes from 
scales in the linear regime, where PT is believed to be 
valid, making the corrections from the nonlinear scales 
small, although probably not negligible. We will try to 
quantify this in more detail below. We will argue that 
halo model requires Rh to be a typical halo radius. 

The philosophy we will advocate in this paper is that 
any analytic approach must give reliable results both in 
Fourier and in configuration space. Failure to do so is 
a sign of something missing in the model. For exam¬ 
ple, a Taylor series may give reliable results in Fourier 
space up to a certain fc, but if truncated at a certain or¬ 
der it generally diverges at high k and makes its Fourier 
transform imposible to calculate. The reason is that the 
series is not convergent at high k, and one has to adopt a 
different summation of the terms that have a better con¬ 
vergence. Our goal is to develop a model that is rooted 
in PT as much as possible, but is also able to reproduce 
simulations. Since simulations have been verified at 1%, 
we will strive for this precision in this paper. 


II. HALO ZELDOVICH MODEL AND 
PERTURBATION THEORY 

The main ingredients of our approach are the following: 

1) 1 loop SPT has loop integrals which, for low k and 
high z, are entirely in the linear regime and thus reliably 
computed. There is no guarantee that the entire loop in¬ 
tegral is correct at all redshifts, even for low k, but we will 
make this assumption here and derive the consequences. 
We will therefore assume that EFT corrections to 1 loop 
SPT are negligible at low k. For higher k and low z SPT 
predictions become increasingly unreliable and will not 
be used. Similarly, 2 loop SPT integrals are negligible 
at high z, while for low z they extend deeply into the 
nonlinear regime and are grossly overestimated in SPT 
[5]. Here we will simply ignore 2-loop SPT, with one 
exception, discussed next. 

2) Zeldovich approximation gives approximately cor¬ 
rect physical picture of how the particles are displaced 
up to the process of halo formation, which stops the 
particles from displacing. The latter has very little ef¬ 


fect on the Zeldovich displacement: most of the displace¬ 
ment is generated by modes in the linear regime and we 
will not be correcting Zeldovich approximation. In terms 
of SPT Zeldovich approximation receives contributions 
from loops at all order, but only from very specific terms 
related to the linear displacement field correlation func¬ 
tion. 

3) Halo formation has to be an essential part of the 
complete model. Halos are objects of very high density, 
leading to a nearly white noise like contribution to the 
power spectrum at low k, with the halo profile parame¬ 
ters determining deviations from white noise at higher k. 
At high k, the halo term contribution dominates the cor¬ 
relations: all of the close pairs are inside the same halos. 
At a scale k ~ l/2i?vir the number of close pairs involves 
all of the pairs inside the virial radius, which must give 
a contribution of the other of M^, where M is the virial 
mass of the halo. Integrating over all the halos, and 
weighting by the halo mass function dn/dM, one obtains 
an estimate of one halo amplitude p~^ J M‘^(dn/dM)dM, 
where p is the mean density of the universe. At redshift 0 
the integral is dominated by cluster mass halos with virial 
radius of l-2h“^Mpc. We expect the one halo term to be 
approximately of this amplitude at A: ~ 0.2 — 0.4h/Mpc. 
However, the mass has to be conserved so the halos have 
to be compensated, which forces the one halo term to 
vanish at very low k. There is no unique way to do this, 
since it depends on what we compensate against. The 
compensation is by definition a two halo term, since the 
mass is being compensated by the particles outside the 
virial radius. Here we compensate against Zeldovich and 
demand that the total agrees with SPT, which then au¬ 
tomatically enforces mass and momentum conservation. 

4) More specifically, we will assume that 1) can be con¬ 
nected to 2) and 3) at some low k, which is low enough 
that SPT can be assumed to be valid, yet large enough 
to still be close to the scale which dominates the com¬ 
pensation. We will match the two on this scale where 
both descriptions are valid, and use 2)4-3) at higher k. 
Since we believe that Zeldovich approximation is a good 
starting point for any modeling we will include it as one 
ingredient of the theory, and decompose the power spec¬ 
trum and correlation function into two parts, 

P(fc) = Pzci + Pbb, Cir) = izeiir) + (1) 

Here subscript Zel stands for Zeldovich and BB for 
broadband beyond Zeldovich [S], which is our one halo 
term. While there may be a residual BAO wiggle sig¬ 
nature that is not captured by Zeldovich, it is essen¬ 
tially negligible in the power spectrum and at most a 
few percent in the correlation function around BAO 
(lOOMpc/h), probably too small to be observed by ex¬ 
isting or future redshift surveys due to large sampling 
variance errors. Here we will thus focus on the model¬ 
ing of the broadband one halo component and ignore the 
wiggle part. Note that we do not assume that the Zel¬ 
dovich part is uncorrelated with the one halo term. In 
this sense our one halo term is not the so called stochastic 
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term uncorrelated with Zeldovich. 

As we argued the key physics ingredient missing in PT 
is the halo formation, which leads to a large contribution 
from the near zero lag correlations, also called the one 
halo term in the halo model [Mia. The halo model 
postulates that the nonlinear evolution leads to halo for¬ 
mation, and that all the dark matter particles belong to 
collapsed halos, with the halo mass distribution given by 
the halo mass function dn{M). The correlations between 
the dark matter particles can be simply split into cor¬ 
relations within the same halo, the one halo term, and 
between halos, the two halo term. On large scales the 
latter reduces to linear theory The one halo term 

in the halo model has a simple physical interpretation on 
small scales, which is that all dark matter particles are 
inside the dark matter halos distributed with a radial 
halo density profile, and this leads to a power spectrum 
that is simply an integral over the halo mass function 
times the Fourier transform of the halo profile squared. 
The halo profile has a compact support, extending out 
to roughly the virial radius, and its correlation function 
is a convolution of the profile with itself, extending to 
roughly twice that. In power spectrum the convolution 
becomes a square of the Fourier transform of the profile, 
which can be expanded as a series of even powers of k 

m, 

PBBik) = Fik)Ao (1 - + •■•) ■ ( 2 ) 

Here the parameters AgRih^n have a specific interpreta¬ 
tion in terms of the integrals over the halo mass function 
in{M) times halo mass M squared, and times 2n mo¬ 
ments of the halo radius averaged over the halo density 
profile [H]. Specifically, Aq = p~^ J M‘^dn{M) is just a 
weighted halo mass squared divided by the mean den¬ 
sity p and does not depend on the halo density profile. 
These arguments however do not yet account for the halo 
compensation. Above we introduced F{k), which is the 
compensation function, required to vanish in fc —>■ 0 limit 
as long as the two halo term converges to linear theory 
in the same limit: mass conservation requires that the 
leading nonlinear one halo term cannot be a constant Aq 
[Hj . Thus one halo term has to be generalized to include 
mass compensation effects: nonlinear effects cause the 
dark matter to collapse into dark matter halos, bringing 
in mass from larger scales, so it has to be compensated by 
a mass deficit at large scales to satisfy the mass (and mo¬ 
mentum) conservation. Because of this one can show that 
the one halo term has to scale as at low k [iniiiiiiii]. 

The two halo term can also be expressed as a convo¬ 
lution of the linear theory over the halo profiles, and the 
resulting Taylor expansion is given by a similar series 


leading order correction the usual EFT term |5]. It is 
clear that both one halo and two halo Taylor expansions 
break down for k > R~^\. The breakdown of the two halo 
term does not matter: at the relevant k the correlations 
are dominated by the one halo term. For the latter how¬ 
ever, a different expansion is advantageous, as we discuss 
below. 

In this paper we argue that the natural way to con¬ 
nect SPT to small scale nonlinear effects is in the context 
of the Zeldovich approximation plus a compensated one 
halo term. In this picture we can think of Pzei as the 
leading order two halo term, and Pbb as the one halo 
term. The motivation for this is that Zeldovich correla¬ 
tion function is almost exact for r > 5Mpc/h, and that 
the correction relative to it is negative, suggesting a com¬ 
pensation of a nonlinear term ^bb{i^)- All the corrections 
to the Zeldovich model thus go into the compensation 
term F{k). In the halo model these corrections would 
arise from compensation of the halo term and from two 
halo term correlations of particles inside the halos. We 
do not try to separate these into the latter that one ex¬ 
pects to scale as k’^P^ik) at low fc, and the former that 
scales as at very low k. In general it is difficult to 
do this separation, as both of these arise from the two 
halo term correlations. It is also not clear that in the 
regime where it matters (fc > OAh/Mpc) these low k ex¬ 
pansions still apply. We will however use them at very 
low k ^ 0.02h/Mpc, where we expect SPT to be valid. 
However, it should be clear that the term Pbb is not 
just the one halo term in the traditional sense, because 
of these compensation corrections at low k. 

In [H] this function F{k) was simply fitted using a 
polynomial form. We will begin by keeping this func¬ 
tion completely general, and then choosing a very simple 
form for it. Our one halo term, and thus the compen¬ 
sation form of F{k), is defined relative to the Zeldovich 
term. If we adopted a different form for the two halo term 
we would obtain a different form of F(k). In principle, 
Zeldovich approximation itself could contain some halos 
with correct mass and so a part of what we usually call 
one halo term could already be contained there. How¬ 
ever, we will see that this must be a minor effect and the 
one halo term we derive agrees with the expected value 
from the halo model. 

While the one halo expansion in even powers of k above 
works in Fourier space up to the virial radius scale of or¬ 
der k ^ Rih 2 Ih/Mpc, it breaks down above that. 
Moreover, powers of k diverge at high k and do not have 
a well defined Fourier transform, making this form un¬ 
suitable for correlation function predictions without re¬ 
summing it first. Instead we will use in this paper the 
Fade series ansatz 


P2h{k) = Pdk) (1 - k^Rlh,2 + k^PthA-) ■ (3) 

The leading order correction scales as /c^PL(fc).R 2 h 
where R 2 h ,2 is also related to an average second moment 
of halo density profile, although with a different mass 
and halo bias weighting m- Note that this gives at the 


Pbb = 


AoF{k) 


1 + \kR^)^^ 

1+EnTlikRnh)^^ 


(4) 


By requiring the series in the denominator to run to a 
higher power than in the enumerator we guarantee that 
the series does not diverge at high k and has a finite 
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Fourier transform. Here we will explore the truncation 
of the series at rimax = 0,1, 2. In terms of the halo model 
Aq has the same interpretation as before, it is the halo 
mass mass squared averaged over the halo mass function. 
It is the only quantity that has units of power spectrum, 
all the other parameters have units of length. For the 
Rm and Rnh parameters we expect that they will be 
related to a typical scale of the halos. There are sev¬ 
eral halo scales one can define: one is the scale radius 
i?s, dehned as the scale where the slope of the density 
profile is -2. Another scale we can define is the virial 
radius R^w = cRg, where c is the concentration param¬ 
eter with values around 3-4 for the most massive halos 
and increasing towards less massive halos m- The mean 
overdensity at i?vir is 200 by definition, while the typical 
mean overdensity at Rs can be of order 1000 or more. 
In the halo model approach we integrate over these with 
the halo mass function, and the number of pairs for each 
is proportional to , which gives most of the weight to 
the very massive halos. As a result, if we use Umax = 1 
and only have one scale parameter, which we denote Rfh, 
then we hnd that the typical scale, when averaged over 
all halos, is of order IMpc/h at z = 0, a typical scale of 
a cluster. 

While the halo model has been successful as a phe¬ 
nomenological model, its connection and consistency 
with PT has not been explored. In this paper we propose 
a Halo Zeldovich model applied to Perturbation Theory 
(HZPT) approach, in which we connect the halo model 
to PT in the regime where both can be expected to be 
approximately valid. We take the approach that one loop 
SPT has a regime of validity on very large scales and gives 
us the correct description of the onset of nonlinearity. 
This is not guaranteed by SPT: the one loop SPT may re¬ 
ceive contributions from small scales which are nonlinear 
and thus not reliably computed. For CDM type power 
spectra, the integrals are convergent and for sufficiently 
high redshift all of the one loop integral contributions 
come from linear scales, the prediction is reliable, and no 
EFT correction is needed. For now we will simply assume 
there are no corrections to SPT at low k, and return to 
this discussion later. 

Let us therefore assume that one loop SPT is cor¬ 
rect at very low k, and that it can be matched to the 
halo model ansatz in equation in the regime of its va¬ 
lidity. The low k limit of Zeldovich approximation is 
0 Pzci = (1 - + fc'‘o-i/ 2 )PL + Q 3 / 2 , where Pl 

is the linear power spectrum, = I/(67r^) f dqP{q) 
is the square of linear displacement field dispersion and 
Qs = l/(107r^)A:^ J dqP‘^{q)/q^ = C^k'^ is a mode cou¬ 
pling integral, as defined in 0. It has been shown in 0 
that this expansion is valid for k < O.lh/Mpc. Note that 
we kept terms beyond 1-loop in the Zeldovich. One loop 
SPT can be written as Pspt = Pl + P13 + P22 , which at 
low k is Pi3 = —61/105fc^cr|PL and P 22 = dSQs/QS. 

Let us begin by first dropping the 2-loop term 
k‘^a^PL /2 from Zeldovich approximation. Then all of 
the one loop terms scale as the square of the power spec¬ 


trum. Matching at low k gives 
PsPT(fc)-Pzei(fc) = = AoP(fc), (5) 


where the linear order cancels in the difference and we 
have dropped all higher order terms since they are neg¬ 
ligible at low k. What does this imply for the amplitude 
dependence of Aq and F{k)l We will adopt the stan¬ 
dard CTg normalization for the amplitude of fluctuations, 
where a's(z) is the rms fluctuation of spheres of radius of 
8Mpc/h, and which is redshift dependent. Often this 
is phrased in terms of the redshift dependence of the 
growth factor D{z), and in linear theory we would write 
CT8 (z) = D{z)as{z = 0). But the amplitude dependence 
is more general than the redshift dependence, since it en- 
compases the idea that changing the redshift or changing 
the amplitude should give the same result in the context 
of PT. Since both SPT and Zeldovich at low k scale as 
a square of the power spectrum, which itself scales as 
CTg, requiring equation!^ to be valid over a broad range 
of k where F{k) is rapidly changing, there can only be 
one solution to the amplitude dependence, Aq oc ct| and 
F{k) = const. This simple result is in very close agree¬ 
ment with simulations spanning a wide range of redshifts 
and models [T3], where the slope of 3.9 was derived. We 
will see below that including the full Zeldovich instead 
of its lowest order further improves the agreement on the 
slope. Note that this is valid for a general compensation 
function F{k). 

To proceed we need to assume a specific functional 
form for the compensation term. A simple way to achieve 
compensation is to use F{k) = 1 — 1/(1 -|- k^R^), which 
has an analytic Fourier transform and vanishes as /c —>■ 0, 
so that the overall model for one halo term is 


Pbb = Aq 



1 


1 -h 


1 + \kRm)^^ 

l + EnTlikRnh)^'-^ 


( 6 ) 

The parameter R governs the transition of the one halo 
term to 0 at low k: this is the compensation scale pa¬ 
rameter, and we expect it to be very large compared to 
the typical halo size. Here we simply chose the simplest 
form for this function Fik) = (1 — 1/(1 -|- k^R^), but we 
do not expect this to be the correct form at all fc, since, 
as seen from equation 0 the residual between SPT and 
Zeldovich at low k is not but k'^Pi,. We evaluate SPT 
and Zeldovich in the low k limit and fit for parameters 
Aq and R. We hnd there is a range of k where the ht is 
good (hgure[^ and we use 0.02h/Mpc < k < 0.05h/Mpc 
for the hts. The ht is not perfect, a consequence of 
the simplihed form of F{k) term, but is a good ht for 
0.02h/Mpc < k < 0.05h/Mpc. Even for the more general 
forms of F{k) one can still dehne a typical compensation 
scale on which F{k) goes from unity to zero. Current 
form of compensation is sufficient for our purposes and 
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FIG. 1: HZPT model for rimax = 0,1,2 (dotted, dashed and solid lines respectively), together with SPT (dot-dashed line) are 
shown as a function of k. We match SPT to HZPT around k ~ 0.02h/Mpc to derive and R, while the remaining parameters 
of our model {Ri, Rih, R 2 h) are fitted to simulations (points) at scales k > O.lh/Mpc. Gray band shows the change in the 
HZPT, nmax = 0 model in case where R value varies for ±3Mpc/h. This plot is for z = 0, and higher 2 versions are similar. 
Simulation points are taken from [^, and are for ag = 0.807 Lambda GDM model. 


we will not explore more general forms. We find 
^0 = 750 (h/Mpc)^ 

R = 26.0 (Mpc/h) (7) 

This is a remarkably simple result. Moreover, it 
agrees well with the recent numerical determination of 
Aq amplitude from a suite of 38 emulator simulations 
at different redshifts in where a scaling Aq = 

670(cr8/0.8)^-®(Mpc/h)^ has been derived over the red- 
shift range 0-1. The amplitude is a bit different because 
the form of compensation used in this letter is a bit dif¬ 
ferent than in mi, but they both provide equally good 
fit to the simulations, as shown in figure The two 
parameters are correlated, and parameter R is less well 
determined than Aq: R comes with five times larger rela¬ 
tive error than the error on Aq. Figureshows the band 
over which R is varied by 12%: a reasonable estimate 
of its error is 3-5% (and the error on Aq is 1% or less). 
We find the slopes of Aq and R to be uncertain at 0.1 
level. The value of the amplitude also agrees very well 
with the expected amplitude of the one halo term in the 
halo model, which is p~^ f M'^dn, where M is the halo 
mass, dn/dM is the halo mass function and p is the mean 
density of the universe. This suggests that Zeldovich ap¬ 
proximation on itself does not contribute much to the one 
halo term. 

To qualitatively derive a value for R and let 
us look at the low k limit of SPT, focusing on 
the leading order and k"^ terms at the peak 
of the power spectrum around fcpeak ~ 0.02h/Mpc, 
where T’L(fcpeak) ~ 27000(cr8/0.8)^(h/Mpc)^ and cr| ^ 


36((T8/0.8)^(Mpc/h)^. Since the difference between Zel¬ 
dovich and SPT for the low k gives (44/105)/c^cr|PL, by 
equating that to k'^AQR? (low k limit of equation]^ we 
derive R ^ (0.42cr|PL(fcpeak)Mo)^^^ 24Mpc/h for the 
best fit value of Aq, in a qualitatively good agreement 
with the value of 26Mpc/h derived numerically. Thus 
the value of R is determined by the linear power spec¬ 
trum amplitude at the peak, rms displacement field and 
the amplitude of the one halo term Aq. To determine 
both Aq and R we need to expand SPT and Zeldovich 
up to k'^ around PL{kpea,v)^ for which we also need to 
numerically evaluate Ca ^ 2 x 10®Mpc/h. Matching 
terms gives [(1/2 - 45/98)^3 + {alfPL/2] = AqR'^, 
which, when combined with k^ term gives R ~ 20Mpc/h 
and Aq ^ 700(Mpc/h)^ at z = 0. We note further that 
there is a considerable contribution beyond from Pl3 
already around k ~ 0.02h/Mpc, and in the numerical fits 
there is some correlation between Aq and R. 

The Zeldovich term beyond 1-loop, [<7\)'^Pl/2, cannot 
be neglected compared to the rest of terms even for 
k < 0.05h/Mpc, and as a consequence the scaling of Aq oc 
ct| and R = const is mildly broken. This term is larger for 
low redshifts, explaining why the slope of ag, scaling is less 
than 4.0 and closer to 3.8. Remarkably, that is exactly 
what simulations suggest. Given the uncertainties in the 
form of F{k) we cannot address in detail the remarkably 
small scatter of Aq against the amplitude an,3 when the 
shape of the power spectrum is varied El- 

The leading order for our compensation term of the 
one halo F{k) is It is often stated that the mass and 
momentum conservation effects generate k^ tail [18] , and 
so one would naively expect the leading term of one halo 
to be k^. This term is generated by P 22 in SPT. How- 
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ever, nonlinear evolution also leads to propagator effects 
contained in ^3, which at the leading order give k^PL 
corrections to the linear theory and we have seen that 
these effects dominate at low k. We have assumed the 
two halo term to be Zeldovich approximation, which con¬ 
tains part of the term in SPT but not all, and so 
the difference is still given by k^Pj^ term at the leading 
order, which has to be the leading order one halo term, 
and happens to coincide with k^ at the peak of the power 
spectrum around k ^ 0.02h/Mpc, where we fit to our 
ansatz. As discussed above, we could use a more general 
form of F{k) that would make the compensation and two 
halo terms exact, but we found this makes no practical 
difference to the final results. It is also possible to make 
a different two halo ansatz where k^PL is entirely can¬ 
celled, for example, our two halo ansatz could be simply 
Pl3 of SPT, or its nonlinear propagator version [18]. As 
shown in [18] . this ansatz leads to one halo term that is 
several times larger than expected in the halo model. It 
is possible to generate a valid halo model based on a dif¬ 
ferent ansatz for the two halo term, but we do not pursue 
this further here. 


III. POWER SPECTRUM PREDICTIONS OF 
HZPT 

In this section we ht higher order parameters of the one 
halo term to obtain the best possible agreement against 
simulations. The best fits for these parameters as powers 
of amplitude crs{z) give 

< = l.87Mpc/h j 

/ ( W 0-29 

i?i„ = 3.87Mpc/h j , 

/ xO.88 

Ri = 3.33Mpc/h^ , 

R 2 h = I.69Mpc/h j . (8) 

Here refers to n„iax = 1 case, while Rit, R 2 h and 
Ri refer to n„iax = 2 case. Even though we only fit to 
one set of simulations, we expect that these parameters 
are nearly universal and apply well to all cosmological 
models, just as in the case of the halo plus Zeldovich 
model of [14]. The main difference relative to [14] is the 
form of the compensation function T'(fc), which was fitted 
to a 10-th order polynomial in m, while here we adopt 
a much simpler form of equation [6 and the expansion of 
the one halo term, which was in even powers of k in |14| , 
while we use Pade expansion here. 

In figure we show the HZPT for P{k) against the 
N-body simulations for four different redshifts. We find 
that the model with Umax = 1 can fit the simulations at 
1% up to fc ^ 0.3h/Mpc, while for rimax = 2 the fits are 
good to 1% up to fc ^ Ih/Mpc. Parameter R ensures 


the low k behaviour of the model, while Aq sets the peak 
amplitude of Pbb which is around k ~ 0.12h/Mpc for all 
redshifts. In [TH] has been argued that SPT is a relatively 
good description against N-body simulations for z > 3. 
This is consistent with our results: as shown in figure]^ 
for z = 2 SPT differs from simulations by only 2%, and 
this difference is presumably even smaller for higher z. 

We have established that one loop SPT can determine 
Aq and R. Which other coefficients of equation can 
PT determine? Since coefficients Rnh, Rn, n > 0, are 
determined by a typical halo scale averaged over the 
halo profile and averaged over the halo mass function, 
they depend on the regime where the overdensity is very 
large: halo virial radius i?vir is defined at a mean over¬ 
density of 200, and Rg is even smaller (with the corre¬ 
spondingly larger mean overdensity). These cannot be 
computed from standard PT, which is only supposed to 
work in the regime where <5 < I. One can show (fig¬ 
ure that these terms become important at a percent 
level around k ~ 0.2h/Mpc at z = 0. It seems un¬ 
likely that PT can make predictions at a percent level 
for k > 0.2h/Mpc, regardless of which PT formalism we 
use. It is however remarkable that one loop SPT can 
predict the amplitude of the one halo term Aq. This is 
possible because Aq depends on the total halo mass that 
has nonlinearly collapsed and not on its profile. At low 
k the non-perturbative halo profile effects also give rise 
to the two halo correction term of order where 

R 2 h ^ IMpc/h. This is the EFT term of [8]. At low k 
this term is a small, a few percent, correction to the SPT 
term, which is of course small compared to linear theory. 
At higher k this term is modified and in our model ab¬ 
sorbed into the overall compensated one halo term Pbb ■ 


IV. CORRELATION FUNCTION 
PREDICTIONS OF HZPT 

We would like to require from a good PT model that 
it works both in Fourier and in configuration space, but 
so far there has been no successful model achieving this. 
Typically, Zeldovich approximation works quite well in 
correlation function but fails in the power spectrum, 
while SPT does not give very good correlation function 
predictions, specially around BAO. In HZPT approach, 
we expect Zeldovich term to dominate the correlation 
function at large radii. In the absence of compensation 
the one halo term would be limited to scales around twice 
the virial radius and below. With compensation these ef¬ 
fects extend to large radii, but as we will show, remain 
small. There is thus a crucial difference of the effect of 
the one halo term between correlation function and power 
spectrum: the one halo term is mostly a few percent ef¬ 
fect for r > 5Mpc/h in correlation function, caused by 
compensation effects. On the other hand, one halo term 
can be very large for k > 1/R ~ 0.04h/Mpc in the power 
spectrum and dominates P{k) for k > 0.2h/Mpc. 

An advantage of the ansatz in equation is the exis- 
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FIG. 2: Simulation power spectrum (points) is shown relative to HZPT model prediction (Pmodei) with rimax = 2 (black solid 
line) at four different redshifts {z = 0.0, 0.5, 1.0 and 2.0). Also shown are SPT (dot-dashed line), HZPT with Umax = 1 (dashed 
line) and linear theory (dotted line). Simulation points are taken from [^, and are for ag = 0.807 Lambda CDM model. 


fence of analytical Fourier transform for low rimax- On 
scales r ^ i?haio, and assuming R ^ Rnh,Rn, one can 
use rimax = 0 and find 




Ape R 


(9) 


Similarly, if we keep the leading Rhi effects in rimax = 1 
case we get: 


iBB{r) 


Ape-^ / ( RV \ R + Rhi' 


( 10 ) 

First equation effectively reduces the low k model that 
we started with in Fourier space to a model in the conhg- 
uration space with the same two parameters, Ap and R, 
that have been determined using SPT. Since the virial 
radius of the largest halos is about 2Mpc/h at z = 0, 
we expect the transition between the two regimes to be 
around 4Mpc/h. The results for the correlation function 
are shown in figure We see that our model signifi¬ 
cantly improves upon other PT results, achieving 1-2% 
agreement down to 5Mpc/h at all redshifts. A possible 
exception is the BAO wiggle, r>80Mpc/h, where there 
may be additional wiggle contribution that was discussed 
in [^, but it is unclear whether it is real given the large 
sampling variance fluctuations. It also improves upon 
Zeldovich approximation, which is already very good by 
itself, with only a few percent deviations from simula¬ 
tions over this range. Our model reduces the correlation 
function relative to Zeldovich, as expected by the com¬ 
pensation effects, which take mass from large scales to 
enhance one halo terms on small scales. As expected our 
resuts agree with SPT on large scales, but only away from 


BAO, making the range where SPT agrees with simula¬ 
tions and our model at 1% level only around 50-70Mpc/h. 

Given that the one halo contribution relative to Zel¬ 
dovich is a few percent only, any further corrections have 
a very small effect on ^(r) for r > 5Mpc/h. For example, 
one can improve the agreement of the model with simu¬ 
lations somewhat by increasing R. We can do this at a 
few percent level since this is the formal error from the 
fits to SPT. In the context of our approach, increasing 
R beyond its SPT predicted value at a level more than 
this would only be possible if there are nonperturbative 
or 2-loop corrections at low k, since an assumption of 
our model is that one loop SPT is the correct theory at 
low k. Note that EFT corrections as advocated by [8] 
will reduce ApR?, which goes in the opposite direction. 
Nevertheless, as discussed above, there is no guarantee 
for SPT to be true even for very low fc: the loop integrals 
can extend into the regime where density perturbation 
exceeds unity. At z = 0, for CDM models, this happens 
around fcni ~ 0.1 — 0.2h/Mpc. The low k expansion of 
SPT is dominated by k^a\P\in{k) terms, which depend 
only on the convergence of the a\ integrals. This con¬ 
verges to about 90% of the value for g < fcni ~ 0.2h/Mpc 
at z = 0, and converges even more to its full value for 
z > 0 where k^\ > 0.2h/Mpc [5]. So we expect the one 
loop SPT to be almost perfectly valid at low k for high 
redshifts, but there may be low redshift corrections at 
a several percent level, which may allow for additional 
changes in R or Ap beyond the predicted value at at 
comparable level. In addition, our form of the compen¬ 
sation term is just an assumed ansatz, which could be 
modified for a better agreement. 

A related question is whether one should include higher 
loop contributions. At a next order in PT is the two 
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FIG. 3: Simulation correlation function (points) is shown relative to HZPT model prediction (^modei) with n^ax = 0 (black 
solid line) at four different redshifts {z = 0.0, 0.5, 1.0 and 2.0). Note that HZPT ^modei with rimax = 0 if completely determined 
with the Ao and R parameters obtained from the low k of SPT power spectrum. Also shown are HZPT with Umax = 1 (dashed 
line), Zeldovich (long dashed line), SPT (dot-dashed line), and linear theory (dotted line). Simulation points are taken from 
[5], and are for ag = 0.807 Lambda CDM model. 


loop SPT or, similarly (although not equivalently), one 
loop LPT. In both cases the corrections to one loop SPT 
can become important at low redshifts, even at low k. 
For example, at z = 0 the two loop SPT correction to 
Pi 3 is about 15% [5], and since the relative correction of 
two loop to one loop scales as cr|, the correction is much 
smaller at higher redshifts. However, higher loop contri¬ 
butions are likely to be grossly overestimated in PT. For 
example, comparison against simulations suggests that 
the one loop LPT contribution to rms displacement is 
in reality almost entirely suppressed, such that the total 
nonlinear value of differs by only 1-2% relative to 
the linear value (j\ at all redshifts [anniiii]. Physically 
this can be understood by the process of halo formation, 
which stops particles from displacing on small scales, and 
instead traps them inside the dark matter halos: the large 
scale displacements, which are one loop in SPT, are cor¬ 
rectly predicted, while the small scale displacements, two 
loop and higher in SPT sense, are strongly suppressed. 
It thus seems better to drop two loop terms entirely, al¬ 
though we have no formal proof of this statement. 

In summary, formally one cannot exclude corrections 
at low k, which will be of the EFT form k^P^, but these 
are likely to be of order a few percent only. We see no 
need for such corrections in our approach: our model is 
accurate at the current precision of simulations. We thus 
argue that one loop SPT is close to the correct theory 
for k < R~^, but this is not a result that can be formally 
derived. 


V. CONCLUSIONS 

In this paper we develop a model for dark matter 
power spectrum and correlation function that is 1% ac¬ 


curate for both, and that is based on perturbation theory 
(PT) as much as possible. We argue that PT approaches 
to large scale structure can only have a hope of being 
valid for very large scales, k < 0.05h/Mpc, a regime 
that we usually do not focus on when comparing PT to 
simulations, since deviations from linear theoryare very 
small there. We also argue that Zeldovich approxima¬ 
tion is a useful starting point for any halo based model, 
and that halo formation has to be essential part of the 
model. We propose a model which we call Halo Zeldovich 
PT (HZPT) model, in which Zeldovich approximation is 
supplemented with the one halo term, and the sum of 
the two is connected to one loop standard PT at low k. 
Within this model we derive the one halo term amplitude 
Aq = 750(Mpc/h)^cr|'®, which agrees with simulations 
both in amplitude and in erg scaling. 

The one halo term needs to be compensated by the 
other halos for the mass conservation, and there are non¬ 
linear contributions from two halo correlations, both of 
which we model this using a very simple functional form. 
This compensation scale of the one halo term has effects 
on the power spectrum at a percent level or smaller, but 
we have argued it is essential in order to have a self- 
consistent model that connects the halo model to PT: its 
introduction gave us one percent accuracy on both the 
power spectrum and the correlation function. In partic¬ 
ular, the deviations of the correlation function of simula¬ 
tions from Zeldovich is negative and a few percent only, 
and this term explains its origin. We have argued that 
the regime where SPT is valid in the power spectrum, 
is at best limited to A: ~ 1/i? ^ 0.04h/Mpc, while for 
higher k one halo term (generalized by the two halo term 
corrections encoded in F{k)) begin to dominate. It is not 
possible to formally exclude presence of nonperturbative 
or higher loop correction terms even at very low fc, but 
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we see no need to consider them and they are likely at 
most several percent for our universe. These terms will 
also generate a correction to the BAO wiggles [5] that we 
have ignored in this paper. 

We have proposed a Bade type expansion of the one 
halo term as a useful functional form that allows one to 
go beyond the convergence radius of the Taylor expan¬ 
sion. We have argued that the value of this radius is 
around IMpc/h, a typical virial radius of halos properly 
averaged over the halo mass function, hence Bade ex¬ 
pansion is necessary if one wants a valid description for 
k > Ih/Mpc. Our approach is similar to the treatment 
of nonlinear redshift space distortions (the so called Fin¬ 
gers of God, FoG), which also require to have an expres¬ 
sion valid for k > 0.2h/Mpc (the FoG scale is typically 
5Mpc/h). In the context of FoG a Lorentzian distribu¬ 
tion is often used, which is just the Fade series at the first 
order. Fade expansion also has the advantage of being 
convergent in both the power spectrum and the correla¬ 
tion function, making the Fourier transforms calculable. 
With this expansion, and keeping terms up to 2nd order, 
we are able to match power spectrum to better than 1% 
against simulations, up to fc = Ih/Mpc. The correlation 
functions also agree against simulations to this accuracy 
down to 5Mpc/h. We expect that a Fade ansatz for one 
halo term will be useful in modeling other correlation 
functions as well, such as galaxy-dark matter and galaxy- 
galaxy correlations. In principle a Fade expansion would 
also be needed for the two halo terms that scale as Pl, 
but in practice the two halo term is irrelevant on scales 
where this would make any difference (fc ~ Ih/Mpc), so 
a simple Taylor expansion giving rise to fc^P^... terms 
suffices. 

We have argued that the predictive power of PT be¬ 
yond Zeldovich has been reduced to two numbers, Aq and 


R. To improve the model further one needs to provide 
information on the halo profiles in the deeply nonlinear 
regime, which is unlikely to be predictable by the PT. 
It has been argued in [T3] that these coefficients are also 
not predictable by N-body simulations, due to the bary- 
onic effects, so PT is not necessarily inferior to N-body 
simulations. For example, there are baryons inside the 
dark matter halos in the form of gas and stars and these 
can redistribute the matter inside halos beyond what the 
N-body simulations can predict. Baryon gas has pres¬ 
sure and this already changes the total matter profiles 
significantly, and even more dramatic effects arise from 
some feedback models where gas is pushed out of the 
halo center, possibly even dragging dark matter along 
These processes will change parameters associated 
with the halo profile, such as Rnh and i?„, (models of 
p2] suggest these coefficients change at a level of 5-10% 
M). but because of mass conservation Aq changes a lot 
less M- Moreover, most of the cosmological informa¬ 
tion content is already in Zeldovich term and Aq [T3] . So 
while HZPT, and PT in general, may only be able to de¬ 
termine two numbers beyond Zeldovich approximation, 
this may also be all that can be reliably extracted from 
N-body simulations at a two point function level. 
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